DAMTP 95-13 



To appear in Communications in Mathematical Physics 

Tetrahedral and Cubic Monopoles 

Conor J. Houghton 
and 

Paul M. Sutcliffe* 

Department of Applied Mathematics and Theoretical Physics 

University of Cambridge 
Silver St., Cambridge CBS 9EW, England 
c.j.houghton@damtp.cam.ac.uk & p.m.sutchffe@damtp.cam.ac.uk 

March 1995 

Abstract 

Using a numerical implementation of the ADHMN construction, we compute the 
fields and energy densities of a charge three monopole with tetrahedral symmetry 
and a charge four monopole with octahedral symmetry. We then construct a one 
parameter family of spectral curves and Nahm data which represent charge four 
monopoles with tetrahedral symmetry, which includes the monopole with octahedral 
symmetry as a special case. In the moduli space approximation, this family describes 
a novel kind of four monopole scattering and we use our numerical scheme to construct 
the energy density at various times during the motion. 
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1 Introduction 



BPS monopoles are topological solitons in a Yang-Mills-Higgs gauge theory in three 
space dimensions. The equation for static monopoles is integrable, so that a variety of 
techniques are available for studying monopoles and constructing solutions. Monopoles 
of charge one and two are well-understood, with explicit solutions known, but for higher 
charges the situation is not so clear. Despite the integrability of the equation, explicit 
solutions for charge three and above are known only in the axisymmetric case, which 
corresponds to coincident monopoles. Very recently, some progress has been made in this 
area with existence proofs for a charge three monopole with tetrahedral symmetry 
and a charge four monopole with octahedral symmetry. In this paper, we compute these 
monopoles using a numerical implementation of the Atiyah-Drinfeld-Hitchin-Manin-Nahm 
(ADHMN) construction and display their energy densities. 

When time dependence is introduced, the monopole equation of motion is not inte- 
grable. However, analytical progress can still be made, via the moduli space approxima- 
tion 0, from knowledge of the static monopoles. This has been extensively studied for 
the case of charge two monopole scattering 0, but the extension to higher charges has 
proved a less tractable problem. We have made some progress in this area and present 
our results on a particularly symmetric example of charge four monopole scattering. The 
charge four monopole has tetrahedral symmetry throughout the motion, which is the key 
to our construction of the relevant spectral curves and Nahm data. We use our numerical 
scheme to display the energy density at various times. 



2 Monopoles, spectral curves and Nahm data 

In this paper, we study solutions of the Bogomolny equation 

A$ = -\e^JkF,k (2.1) 

for SU(2) BPS monopoles in H^. Here Di = + [Ai, is the covariant derivative with Ai 
the su(2)-va\vLed gauge potential and Fjk the gauge field. $ is the Higgs field, which is an 
su (2 valued scalar field satisfying the boundary condition 

||$|| = 1 - ^ + 0(^) as r^cx) (2.2) 

where r = \x\, = — |tr$^ and is a positive integer, known as the magnetic charge. 
We shall refer to a monopole with magnetic charge /c as a fc-monopole. The energy density, 
of a monopole is given by 

£ = -itr(A$)(A$) - ^tr(F,,F,,). (2.3) 
The energy is the integral of S over all space and is equal to Snk. 
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Equation ( |2.1| ) may be obtained by dimensional reduction of the self-dual Yang-Mills 
equation in M^, for which there is a well-known twistor correspondence; namely that solu- 
tions of the self-duality equations correspond to certain holomorphic vector bundles over 
the standard complex 3-dimensional twistor space. This correspondence may be reduced 
1^, |, ^ to give that monopoles correspond to particular holomorphic vector bundles over a 
mini-twistor space TT, which is a 2-dimensional complex manifold isomorphic to the holo- 
morphic tangent bundle to the Riemann sphere ie Tr=T(DP^. Moreover, the bundle (and 
hence the monopole) is determined by an algebraic curve in TT, called the spectral curve, 
which must satisfy certain reality and non-singularity conditions. 

The space TT is a fibre bundle over (DP^ with each fibre being a copy of (D. Let ( be the 
standard coordinate on the base space and i] the fibre coordinate, then the three spectral 
curves of interest in this paper are 

r(l/6)=^r(l/3)3 , , , 

C(C--1)^0 (2.4) 

+^^^(C^ + 14C^ + 1) = (2.5) 
r;^ +i36aK^7]C{C - 1) + 3k\C^ + 14C^ + 1) = 0. (2.6) 



In [|I| it is proved that (^) is the spectral curve of a 3- monopole with tetrahedral symmetry 



and that ( p.Sj) is the spectral curve of a 4-monopole with octahedral symmetry. It is the 
monopole fields and energy densities which correspond to these two spectral curves that 
we shall compute numerically in Section 3. In Section 4 we shall prove that ( |2.6| ) is the 
spectral curve of a 4-monopole with tetrahedral symmetry for all a G (— 3~^/^-\/2, S~^^'^\/2) 
where 2hi is the real period of the elliptic curve 

y'^ = 4:{x^ - x + Sa^). (2.7) 



If a = then ( p.6| ) becomes ( |2.5| ), so that at this point the 4-monopole has octahedral 
symmetry. 

Although the spectral curve approach to monopoles is a very useful and powerful tech- 
nique, its main drawback is that the non-singularity constraint, which an algebraic curve 
must satisfy to be the spectral curve of a monopole, is rather formidable to check. However, 
there is an alternative approach to the construction of monopoles which nicely complements 
the spectral curve formulation, in the sense that the non- singularity of the monopole is 
automatic. The ADHMN construction |^, Q is an equivalence between fc-monopoles and 
Nahm data (Ti,T2,T3), which are three k x k matrices which depend on a real parameter 
s G [0, 2] and satisfy the following; 

(i) Nahm's equation 

= -eijk[Tj,Tk] (2.8) 
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(ii) Tj(s) is regular for s G (0, 2) and has simple poles at s = and s = 2, 

(iii) the matrix residues of (Ti, T3) at each pole form the irreducible fc-dimensional 
representation of SU(2), 

(iv) T,is) = -Tjis), 

Equation (i) is equivalent to a Lax pair and hence there is an associated algebraic curve, 
which is in fact the spectral curve. Explicitly, the spectral curve may be read off from the 
Nahm data as the equation 

det{r] + (Ti + iTa) - 2mC + (Ti - m)^) = 0. (2.9) 

In Section 3, we review how to obtain the monopole fields from the Nahm data and 
explain our numerical implementation of this procedure. 



3 Numerical ADHMN construction 



Finding the Nahm data effectively solves the nonlinear part of the monopole construc- 
tion but in order to calculate the fields themselves the linear part of the ADHMN con- 
struction must also be implemented 0, |[. Given Nahm data (Ti, T2, T3) for a /c-monopole 
we must solve the ordinary differential equation (ode) 



as 



Ob j ^ j I IrT'j 







(3.1) 



for the complex 2A;- vector v(s), where Ifc denotes the k x k identity matrix, aj are the 
Pauli matrices and x = {xi,X2,X3) is the point in space at which the monopole fields are 
to be calculated. Introducing the inner product 



(vi, V2) = / vlv2 ds 
Jo 



(3.2) 



then the solutions of ( p.l|) which we require are those which are normalizable with respect 
to (|3.2|) . It can be shown that the space of normalizable solutions to (|3.1|) has (complex) 
dimension 2. If vi, V2 is an orthonormal basis for this space then the Higgs field $ is given 
by 

((s-l)vi,Vi) ((S-1)V1,V2) 
((S-1)V2,V1) ((S-1)V2,V2) 

There is a similar expression for the gauge potential but we shall not need this here. 
In some cases this procedure, which goes from Nahm data to the Higgs field, may be 



(3.3) 
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completed analytically to give an explicit closed form for $. However, the Nahm data 
which we consider in this paper is sufficiently complicated that to calculate a closed form 
expression for $ appears not to be a tractable problem. We therefore turn to a numerical 
implementation of the above procedure, which we now describe. 

The first issue we confront in a numerical approach is to calculate numerical values 
for the Nahm data on the interval s G [0,2]. Although we shall have explicit expressions 
for the Nahm data this is still not quite a trivial issue, since the expressions involve the 
Weierstrass elliptic function p and its derivative. However, we can keep the number of 
calculations of Tj(s) to a minimum by noting that if a fixed step ode solver is used to 
integrate then the Nahm data is required at the same s values for every integration 
of (|3.1|) for all initial conditions and x positions. Therefore we compute, once and for all, 
Tj(s) at 2P equidistant points for s G [0, 2] and store these values, which are then used as 
a look-up table when integrating (|3.1|) by a fourth order Runge-Kutta method with fixed 
steplength ds = 2P~^. The values in the look-up table are computed from the closed form 
expressions using mathematica. 

Let ^{I) denote the space of solutions to (|3.1|) which are normalizable for s in the 



interval /. Then we require a basis for the 2-dimensional space fi([0, 2]). The question we 
now address is how to obtain this basis from solutions of the initial value problem (ivp) 
associated with the ode (|3.1| ). Consider the ivp of ( |3.1| ) at the pole s = 0, which has the 
form 

s'-^ = B.. (3.4) 

where Bs is a regular 2k x 2k matrix function of s G [0,2). This is a regular-singular 
problem so that fi([0, 2)) has dimension N, where is the number of positive eigenvalues 
(counted with multiplicity) of Bq. If N was equal to 2 then we could easily compute a basis 
for fi([0,2]) since it would (almost) be given by a basis for f2([0,2)), which can be found 
by integrating ( p.l|) , as described above, with two different initial conditions. However, for 
all the cases considered in this paper we find > 2, so that the problem requires a little 
more work. By symmetry of the Nahm data, if we consider the ivp of ( p.lD at the pole 
s = 2 (with ds < 0) then we have a similar regular-singular problem involving a matrix 
which again has A^ positive eigenvalues. By integrating this ivp we can compute a basis for 
the A^-dimensional space f2((0,2]). The 2-dimensional space we require is the intersection 
of the above two A^- dimensional spaces ie 

fi([o, 2]) = ^]([o, 2)) n ^]((o, 2]). (3.5) 

To find the intersection of these two spaces is a shooting problem but because the ode 
( p. ID is linear this shooting problem can be reduced to linear algebra as follows. Let 



ui(s), U2(s), .., U]\f{s) denote A^ 2fc-vectors which form a basis for f2([0, 2)) and 

UAr+i(s), Ujv+2('5), .., U2n{s) a basis for ^((0, 2]). Explicitly these vectors are computed by 

solving the ivp at s = and s = 2 with A^ different initial conditions each. Define the 
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2k X 2N matrix 



U 



Ui(l) U2(l) ... U2Ar(l) 



(3.6) 



then we need to find a basis for the 2- dimensional kernel of U ie to solve the matrix equation 



Uw = 



(3.7) 



for w = {wi, ...,W2n) 7^ 0. Numerically this is performed by row reduction of the matrix 
U followed by back substitution. Let w*^^) and w*^^) denote two independent solutions to 
(p.7|), then a basis for fi([0, 2]) is given by 



Ms) 



if < s < 1 



-E?f^+iwfu,(.) ifl<s<2 



(3.^ 



for I = 1,2. To summarize, the above procedure consists in integrating (|3.1|) times from 
each end of the interval [0, 2] to the centre and then finding linear combinations of these 
solutions such that these combinations, which start at each end of the interval, match at 
the centre. 

Given vi and V2 we use the Gram-Schmidt orthonormalization algorithm, with inner 
product ( |3.2| ) (and the integral calculated from the P data values via a simple Simpsons 
rule), to obtain two orthonormal vectors Vi,V2- The Higgs field $ is then computed 
according to (^.3[) and to calculate the energy density we make use of the formula 



S = Am 



(3.9) 



where A denotes the laplacian on IR'^. Numerically we use the above scheme to calculate 
11$ IP on a spatial lattice of M x M x M points and approximate the laplacian in ( |3.9D using 
a finite difference method with a 7-point stencil. This completes our numerical ADHMN 
algorithm. 

Although every stage of our algorithm is a relatively inexpensive computing task each 
must be executed many times to build up a detailed picture of the energy density. To 
produce each of the energy density plots appearing later in the paper we used the values 
P = 50 and M = 31, with {xi,X2,X3) G [—5,5] x [—5,5] x [—5,5]. This means that the 
ODE ( |3.1| ) must be solved to the order of 10^ times to produce one energy density plot. 
Implementing our scheme on a workstation gave a run time of approximately 30 minutes 
to produce the data for each plot. 

The Nahm data which correspond to the spectral curves (|2.4|) and ( p.5|) , of a 3-monopole 
with tetrahedral symmetry and a 4-monopole with octahedral symmetry respectively, is 
given in |l| and we shall make use of it now. The method used to obtain these data is 
reviewed in Section 4 when we shall use it to calculate the Nahm data for the spectral 
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curve ( p. 61 ). Not all the Nalim data given in explicitly satisfies conditions (iv) and (v) 
given earlier. However the properties of the associated spectral curves implies that there 
exists a constant k x k matrix in each case such that conjugation of the Nahm data by 
this matrix produces equivalent Nahm data which does satisfy the conditions, and this is 
enough. Conjugation by a matrix is equivalent to a change of basis for the fc-dimensional 
representation of SU(2) formed by the matrix residues of (Ti, T2, T3) at the s = pole. For 
our purposes it is convenient if this is a real representation and so (if necessary) we make 
a transformation to achieve this. In the case k = 3 the Nahm data from [|l| is equivalent 
to 



(3.10) 





" 
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Ti = 








—z 


T2 = 













—z 
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where 



+ 



and p is the Weierstrass function satisfying 



r(i/6)r(i/3) 



= _ 4 



(3.11) 



(3.12) 



where ' denotes differentiation with respect to the argument. 



With this Nahm data equation (|3.1| ) is equivalent to the set of coupled ode's 



Vi + X^Vi + (Xi + iX2)v2 + izv^ - ZVq = 
V2 - XsV2 + (xi - iX2)vi - izvi + zv^ = 
V3 + XsV3 + (xi + ix2)v4 - izvi - izVQ 
V4 - XsV4 + (xi - ix2)v3 + izv2 - izv^ 

V5 + + (Xi + iX2)v6 + ZV2 + izVi = 

ve - x^ve + (xi - ix2)v5 - zvi + izv^ = 






= 
= 







(3.13) 



where "^i = ^ etc. In terms of the notation introduced earlier we find that the matrix 
Bq has eigenvalues {1, 1, 1, 1, —2, —2} so that N = 4. 

Figure 1. displays the output of our algorithm for this case. The plot shows a surface 
of constant energy density S = 0.20. The tetrahedral symmetry of this surface is clearly 
evident and plots for other values of S close to this one are qualitatively similar. For large 
values of S the surface breaks up into four disconnected pieces centered on the vertices of 
a regular tetrahedron. 

We now turn to the 4-monopole with octahedral symmetry. The Nahm data from ^ 
(after a change of basis) is 





V3{2y- 


-lOy 



XI 



X 



V3{2y 


—6y — 2x 






—6y — 2x 


VS{2y - x) 



-lOy 


v^(2y- 




X 
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Figure 1: Tetrahedral 3-monopole; surface of constant energy density 8 = 0.20. 



To 



V^{2y-x) lOy 

-V3{2y-x) -6y-2x 

Qy + 2x ^/3{2y-x) 

-lOy -V3{2y-x) 



T3 = ^ 



-Ay-3x 

12y-x 

-12y + x 

Ay + Sx 



where 



X = ue 



10p'(sweW4/2) ' ^ 10p'(sweW4/2)^ 



and p is the Weierstrass function satisfying 

Then equation ( ^TTl) is equivalent to 

■ui + X3t;i + (xi + ix2)v2 + {4:y + Sx)^;! + 20yv8 = 



(3.14) 



(3.15) 



(3.16) 
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V2 - X3V2 + (xi — ix2)vi — (4y + 3x)v2 + 2VS(—2y + x)v3 = 

V3 + X3V3 + (xi + ix2)v4 + 2VS(—2y + x)v2 + (-12y + x)v3 = 

V4 — X3V4 + {xi — ix2)v3 + (12y - x)v4 + 4(3y + x)v5 = 

■Us + X3V5 + {xi + ix2)vQ + 4(3y + x)v4 + {12y - x)v5 = 

ve - X3Ve + {xi - ix2)v5 + {-I2y + x)v6 + 2y/3{-2y + x)v7 = 

V7 + X3V7 + (xi + ix2)vs + 2VS(—2y + x)vq — (4y + 3x)v7 = 

vg - X3VS + (xi - ix2)v-j + 20yui + {Ay + '3x)vs = 0. 

(3.17) 

We find that the matrix Bq has eigenvalues {|, |, |, |, |, — |, — |, — §}, so that N — 5. 

Figure 2. displays the output of our algorithm in this case. The plot shows a surface of 
constant energy density E = 0.14. Note that for the monopole with octahedral symmetry 
a constant energy density surface could have resembled an octahedron or a cube; clearly 
it is the latter. It is therefore more natural to refer to this monopole not as an octahedral 
monopole but as a cubic monopole. For large values of £ the surface breaks up into eight 
disconnected pieces on the vertices of a cube. 




Figure 2: Cubic 4-monopole; surface of constant energy density £ — 0.14. 



9 



4 Four Monopole Scattering 

In this Section we shall use the method of |l| to construct a one parameter family of 
Nahm data, which represent four monopoles with tetrahedral symmetry. The imposition 
of tetrahedral symmetry facilitates the solving of Nahm's equations, since we will only 
have to consider Nahm data which are invariant under the action of the tetrahedral group 
T C ^0(3). 

The Nahm data are an M'^ (8> sl{k, (D) valued function of s, which transform under the 
rotation group SO (3) as 

3®sl{k) (4.1) 

where r denotes the unique irreducible r dimensional representation of su{2). Since gl{k) = 
k®k_ Clebsh- Gordon decomposition gives 

gl{k) = 2k -I © 2k -3 © ... © 3 © 1 

and so 

sljk) = 2fc-l © 2fc-3 ©...©3 (4.2) 

Substituting into ( |4.1| ) 

3®sl{k) = 3©( 2A;-l © 2A;-3 ©...©3) 

= ( 2fc + 1 © 2A; - 1 © 2A; - 3 ) © . . . © (5 © 3 © 1). (4.3) 

Thus the Nahm data corresponding to four monopoles are in the carrier space (9© 7©5) © 
(7 © 5 © 3) © (5 © 3 © 1). The 1 representation is, of course, invariant under all of 50(3). 
We will calculate this and the other tetrahedral invariants. 

Write X,Y and H for the basis of su{2) satisfying the commutation relations 

[X,Y] = H- [H,X] = 2X- [H,Y] = -2Y. (4.4) 

These may be represented by the principal su{2) subalgebra of sl{k) which in turn acts 
on the algebra by the adjoint action. In this representation X is a rank k — 1 nilpotent 
element and a basis of sl{k) can be generated by acting with Y on X^, for r = 1, 2, .., A; — 1. 
Thus 

X^'-i {adY)X^-^ {adYfX^- 
X"- {adY)X'' {adYfX"" 

X {adY)X {adYfX 

is a basis of sl{k). The element X''" of the abelian nilpotent subalgebra (X, X^, . . . X^^^) 
is the highest weight vector for the su{2) representation 2r + 1 lying in the decomposition 
(13) of sl{k). 



(adYf'-^X''- 



(adYV^'X'- 
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It is convenient to exploit the representations of su{2) on homogeneous polynomials 
over (DP^ since the invariant homogeneous polynomials are known P], and also it connects 
with the spectral curve approach. The r + 1 dimensional su{2) representation r + 1 is 
defined on degree r homogeneous polynomials by the identification 

O r\ 

X = Cij-; Y = Coj-; H = -(oj- + (^j- . (4.5) 

^Co ^Ci ^Co ^Ci 

In the case of degree r homogeneous polynomials we can identify highest weight vector 
C[ and basis (Co^)Cr) • • • ? (Co^)''Cr)}- Thus we can relate a degree 2r homogeneous 
polynomial q2r{Co, Ci) cind a matrix S in the 2r + 1 representation of the decomposition of 
sl{k) by rewriting g'2r(Co5Ci) ?2r(Co^)Ci'' then letting 

S = q2r{adY)X\ (4.6) 

The lowest degree T-invariant polynomial is of degree 6. It is 

Ci'Co - CiCo'- (4.7) 
There is a degree 8 T-invariant polynomial 

Cf + MCfCo + Co' (4.8) 

which is also invariant under the octahedral group. Thus, in addition to the SO (3) invari- 
ant there are T-invariant Nahm triplets lying in the 9 representation and in both the 7 
representations. It is convenient to write (9„ © 7^ © 5^) © (7^ © 5„ © 3^) © (5^ © 3^ © 1^) 
so that we may distinguish 7^ and 7„. 

We can now construct the T-invariant Nahm triplets in 9„ and 7„ by (|4.6| ) and the 
inclusion 



2r+l ^ 3(g) 2r-l = 2r + l © 2r-l © 2r-3 
P2r-(Co,Ci) ^ ii®^^ + ^ioii®j7-T7-+io®^rpr- (4.9) 

We choose matrices 

" V3 

2 

~ V3 



Polarizing ( [4.7D yields 

® (20C?Co) + 26^0 © (5Cf - 5Co') - ^0 ® (20CiCo) (4.11) 



Y 

















' 3 








" 













H = 





1 











2 














-1 


















. 








-3 . 



(4.10) 
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which we put into the form 



5(Co|^)C? + «„|^)??®|5-|j(Co|^)1Cf 



+ 



3/-4 



and convert to k x k matrices 

X ® SarfyX^ + adYX 
These matrices were calculated explicitly using maple and are proportional to 



[5 - ^{adYy]X^ + liadYfX [-l(adY)']X' 
24 2 





" 2V3 

















-2^3 
















Similarly polarizing (^.8|^ 


yields 






e?®(56Cf + 168C?Co) + 


26^0 ^ 






















(224Cf Co') + ^0 





Vs 






(168CfCo+56Co'; 









(4.12) 



(4.13) 









-v^ 
(4.14) 



(4.15) 



which becomes 



d 



1 



7 

15^^°aC: 
d 



d 



28 



2(^0^)^^? 



^^1 

n2 , 



7 



15 
d 



90^^°9Ci 



d_ 

'dCi 

\6i/-6 



'28 d 



(4.16) 



yielding invariant Nahm triplet 












-20 
















-12 














4^3 






Y. 
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Y, 



-2^3 

6 







-4 







-2^3 






12 









-12 

4 



(4.17) 



We calculate the T-invariant in Zm by constructing an isomorphism between it and 7„. 
We observe that X^X^ is a highest weight vector of the representation 2r + 3 ^^ and a basis 
can be generated by successive application of {adY (S) 1 + 1 <^adY). Thus, for example, the 
invariant ( [4. 131) can be written 



[5{adY (g) 1 + 1 (g) adY) - — {adY 



1 + l®adYf]X®X^. 



(4.18) 
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The highest weight vector for 2r + 3 ^ is easily calculated by noting that it is annihilated 
by {adX ® 1 + 1 ® adX). It is {adYX ® X'^+^ - ® adYX'^+^). We can then map 



2r + 3 , 



2r + 3. 



by 



X (8) X'' ^ (adYX ® 
Thus the 7^ invariant is 



-X(8)adyX''+i). 



(4.19) 
(4.20) 



[5{adY ® 1 + 1 <8) ady) - ^{adY ® 1 + 1 ® ady)^](arfrX ® - ® adrX^) 



24 



10 



36 



X (8) {--{adYy + —{adYy)X'^ + adYX (g) { — {adY) + —{adYf)X 
+(adYfX ® (5 - ^{adY)^)X^ 



(4.21) 



with corresponding matrices 






-6 



-2\/3 
6 








-2v^ 







2v^ 




2^3 

-2v^ 



-2^3 












3 


V3 














-3 




















(4.22) 



The easiest way of calculating the 5'0(3) invariant is to observe that it is annihilated 
by both {adX 1 + 1 (g) adX) and {adY 1 + 1 (g) adY). It is 

X ® [adY fx - adYX ® adYX + [adYfX®X = X® (-2F) - adYX ®H + {adYfX(»X 

(4.23) 

We now change basis so that the SO (3) invariant is given by (pi, p2, ps), the su{2) basis 
satisfying [pi, P2] = 2ps etc; 



Thus 



Pi = X-F; p2=tiX + Yy, ps = tH. 



(Fi, Y2, Ys) ^ {Y;, Y^, Y^ = i^Yi + n, ~Yi + lYs, -lY^) 



(4.24) 
(4.25) 



and similarly for (^1,^2,^3) and (VFi, 1^2, W^s)- We drop the primes on the transformed 
quantities. 
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With a view to calculating the Nahm equations the commutation relations satisfied by 
the invariant Nahm vectors were calculated using maple 

[Fi, F2] = -48p3 - Sn, [^1, Z2] = |P3 + i>3 

[W^l, W2] = f P3 + 1^3, [Pl, Y^] + [Yi, P2] = -61^3 

[pi, Z2] + [Zi, P2] = -4^3, [pi, 1^2] + [W^i, P2] = 2W^3 (4.26) 
[Yi, Z2] + [Zi, ^2] = -32Z3, [Zi, 1^2] + [W^i, Z2] = 

Writing 

Ti{s) = x(s)pi + y{s)Yi + + (4.27) 

the Nahm equation ^ = 22] reduces under the requirement of T-symmetry to the set 
of coupled nonlinear equations 

^ = 2x'-A^y' + \z'+^^w\ (4.28) 
ds o o 

^ = _8^2 ^ 3 2 ^ 6 2 _ g^^^ 

as 5 5 

— = -Axz-?,2yz, (4.30) 
as 

— = 2xw + 16w|/. (4.31) 
as 

Calculation of the polynomial det(77 + (Ti + iT2) - 2iT3C + (Ti - iT2)0 gives the spectral 
curve 

r]' + ci77C(C' - 1) + ciC2(C' + 14C' + 1) = (4.32) 



where 
and 



ci = 288^(^2 + Ay^ + 3^^ - 4xy) = 288ic; (4.33) 

48 48 

C2 = (60y2 + 2,z^ - 3w2 + 20xy) = ici (4.34) 

288^^ ^ 288 ^ ^ ^ 

are constants. 

To solve these equations, we observe that w can be set to zero, so we do so. We let 
u — X — 2y and v — x + 8y to get 

^ = 2uv (4.35) 

OS 

;2 = i4 (4-36) 
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and 

, U n cy 3Cl 

Cn = — V-U'' 

Define = — IGc'^^Cg and a = Sc'i/n^ to obtain 



(4.37) 



(111 

4— = -V64m4 - 4/t4 + 3a2KS-2. (4.38) 
ds 



Let t = Ks and m(s) = —kJ p{t)/2 giving 



rfp(t),2_4^4^^^^)2_^^ 3a^)_ ^439) 



rft p(t) 

Thus pit) is the Weierstrass function satisfying 

p'2 =4p3_4p + i2a2. (4.40) 



Hence (^:28|)- (|4:3T|) are solved by 



^ ^(-2y^.i^) ,4.41) 
. ^(v^.ifl) (4.42, 

= 4.43 

2p{Ks) 

w{s) = (4.44) 

In order to determine that these Nahm data correspond to a monopole we need to 
examine the boundary conditions. As t — > 0, p{t) ~ and so 

X ~ — y ~ 0, 2: ~ 0. (4.45) 

Therefore at s = the residue of Tj is — and so it forms an irreducible representation 
of SU(2). As t ^ 2uJi (the real period of the elliptic function p{t)) 

3 1 

X ~ -, y ~ — -, ~ (4.46) 

10s 10s 

where s = 2ujik^^ — s, and so the residue of T, is Ri = —j^Pi + jo^i- The eigenvalues 
of 2i_R3 are {3, 1, —1, —3} demonstrating that the -Rj's are an irreducible representation of 
SU(2). Furthermore the functions x,y and z are analytic for t e (0,2u;i). We set k = ui, 
so that the poles occur at s = and s = 2. This demonstrates the existence of a one 
parameter family of monopoles with spectral curves 

7]^ + i36aK^r]C{C^ - 1) + 3k\C^ + + 1) = 0. (4.47) 
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A single monopole with position (xi,X2,X3) has spectral curve 



r]-{xi+ 1x2) + + {xi - ix2)C'^ = 0. (4.48) 

The product of four spectral curves corresponding to four monopoles positioned at the 
vertices 

{(+6, +6, +b), i+b, -6, -b), i-b, -b, +b), i-b, +b, -b)} (4.49) 
of a regular tetrahedron (where 6 > 0) is 

- 16ib^r]{C' - C) + 46^(C^ + 14C^ + 1) = 0. (4.50) 



The spectral curve ( [4.47| ) has this form when 

a = -3-^/^V2. (4.51) 
Examination of the integral expression for k 

K = r '^'^ (4.52) 

io VI - + 3a2x6 ^ ' 

where X is the first positive real root of = 1 — + Sa^x^, shows that k ^ 00 as 
a ±3-5/^72 but it is finite for a G {-?>~^/^V^,^~^'^V^)- We conclude that as a 
approaches — 3~^/^v^ (|4.47] ) describes the superposition of four well-separated monopoles 
on the the vertices ( |4.49|) of a tetrahedron, with the distance between monopoles equal to 
3^/^^. The tetrahedron dual to the one above has vertices 

{(-6, -6, -6), (-6, +6, +6), (+6, +6, -6), (+6, -6, +b)} (4.53) 

with a corresponding product of spectral curves given by 

T]^ + mb^ri{C^ - C) + 46^(C^ + 14C^ + 1) = 0. (4.54) 

Clearly this is the form of the spectral curve (lOTi) when a = 3-^/^72. 
If a = then z = and k is given by 



dx T{l/Af 



vT- a;' V»7r 



p4 



(4.55) 



so that the spectral curve becomes that of the cubic 4-monopole given by (|2.5|) . 

We have derived a one parameter family of 4-monopoles with tetrahedral symmetry. 
They correspond to a one parameter family of spectral curves. We can use this one param- 
eter family to discuss low energy scattering of 4-monopoles because it is a geodesic in the 
4-monopole moduli space. In order to prove that the one parameter family is a geodesic 
in the 4-monopole moduli space we must allow for the possibility that the function w is 
non-zero. We will find that solutions to the tetrahedral Nahm equations ( [4.28| )-( p:.31| ) cor- 
respond to the same one parameter family of spectral curves irrespective of whether or 
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not w is set to zero. This means that the fixed point set of the tetrahedral symmetry in 
the 4-monopole moduh space is one dimensional and since the fixed point set of a group 
action on the moduh space is totally geodesic this means that this one parameter family 
is a geodesic. 

If w is not set to zero, we find that 

dw 

— = 2wv 
as 



and that the equation for ( |4.35| ) is unchanged. Thus 



d / w 
ds \u 



(4.56) 



which implies that w oc u. For convenience we choose the constant to be ^sinh^. Fur- 
thermore if we set 



then u satisfies exactly the same equations as were formerly satisfied by u and replaces u 
in the expressions for the two constants, c'l and ie. 



du 

— = 2uv 
ds 



V? 



and 



^2 



O '2 

2 ~2 "-""1 
V — U — 



(4.57) 
(4.58) 

(4.59) 



The solutions m, v are identical to the earlier solutions -u, v in the w = case so that 



now 



X 



4m sech 9 + V 



V — u sech 6 



(4.60) 
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w 



■utanh^ 



(4.61) 



It can be seen by explicit calculation that the matrix residues at both ends of the intervals 
are irreducible representations irrespective of the value of Q. Thus the Nahm data always 
corresponds to a monopole. However, it is clear from the above construction that the 
constants c'^ and are independent of Q. Hence changing the value of Q does not change 
the spectral curve. There is a one-to-one correspondence (up to gauge transformations) 
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between monopoles and spectral curves, so changing the value of 6 simply corresponds to a 
gauge transformation of the Nahm data. So, by a suitable choice of gauge we can set = 
without loss of generality. This means that we have arrived at a one parameter family of 
monopoles by forcing the 4-monopole to admit tetrahedral symmetry. This proves that 
the family of monopoles is a geodesic in the 4-monopole moduli space. 

In the moduli space approximation the dynamics of k monopoles is approximated 
by geodesic motion on the fc-monopole moduli space Aik- In this Section we have identi- 
fied a totally geodesic one-dimensional submanifold of Ai^ and so we can use the moduli 
space approximation to convert this into a result on four-monopole scattering. Since our 
submanifold is one-dimensional the explicit form of the metric is not important. The infor- 
mation we lose by not knowing the metric is how physical time is related to the parameter 
a, but this is not serious. The above results, therefore, have the following interpretation 
in terms of monopole scattering. Four monopoles approach from infinity on the vertices 
of a contracting regular tetrahedron, coalesce to form a configuration with instantaneous 
octahedral symmetry, and emerge on the vertices of an expanding tetrahedron dual to the 
incoming one. 




Figure 3: Schematic representation of 4-monopole scattering. 
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To make the above scattering process a little clearer we give a schematic representation 
in Fig 3. We draw a cube whose centre is at the origin and whose edges are parallel to the 
coordinate axes; it is to be associated with the cubic 4-monopole (compare Fig 2.). The 
incoming monopoles are represented by black spheres and the outgoing monopoles by white 
spheres, with an arrow indicating the direction of motion for each. Note that if one tried 
to extend this asymptotic interpretation to the region in which the monopoles are close 
together then one would conclude that the monopoles suffer no deflection and simply pass 
through each other. But this is misleading, since each of the outgoing monopoles cannot be 
identified with a single incoming monopole but is a composition of all the incoming ones. 
(A similar misleading interpretation exists for the scattering of three topological solitons in 
the plane. For k solitons in the plane with cyclic Ck symmetry the solitons scatter through 
an angle ir/k, which for A; = 3 could mistakenly be taken for zero scattering angle). 

To obtain a true picture of the scattering process, one needs to examine the energy den- 
sity during the motion. Using our numerical scheme we can do this. Fig 4.Q shows a surface 
of constant energy density £ = 0.06 for the five values a = —0.25, —0.18, 0.00, +0.18, +0.25. 
We see that, indeed the energy density is initially localized in four regions roughly centered 
on the vertices of a tetrahedron. Let us think of these vertices as being opposite corners of 
a cube as in Fig 3. On any one face of the cube the incoming energy density is concentrated 
on two opposite corners of the face (black spheres in Fig 3.) and it flows around the edges 
of the face until it is localized on the two remaining corners (white spheres in Fig 3.) as 
the monopoles separate. This suggests that a useful way to view this scattering process is 
as pairs of 90° scatterings occuring simultaneously. 

5 Conclusion 

Using a numerical scheme we have computed the energy densities of a 3-monopole with 
tetrahedral symmetry and a cubic 4-monopole with octahedral symmetry whose existence 
was recently proved [l[]. We then proved the existence of a one parameter family of de- 
formations of the cubic 4-monopole which has tetrahedral symmetry. In the moduli space 
approximation this describes a 4-monopole scattering process and we used our numerical 
scheme to analyse this further. 

There are a number of interesting aspects which remain in the study of monopoles with 
the symmetries of regular solids. One obvious task is to construct a family of 3-monopole 
solutions which describes the scattering process in which the tetrahedral 3-monopole is 
formed. This problem is currently under investigation but is more difficult than the scat- 
tering considered in this paper, since the family has only C3 symmetry which is not as 
useful as the tetrahedral symmetry which allowed us to solve the corresponding problem 
in the 4-monopole case. Another issue is that of a monopole with icosahedral symmetry. 
A monopole with icosahedral symmetry has to have charge at least six, but in [|l| it was 

^Fig 4. is not available in the hep-th version of this paper. A hard copy is available on request to 

P.M.SutclifFe@ukc.ac.uk, or it can be viewed at URL 

http://www.ukc.ac.Uk/fMS/maths/people/P.M.Sutcliffe/preprints.htm] 
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proved that no such monopole of charge six exists. We have proved that an icosahedral 
monopole of charge seven exists and are currently investigating its properties. These re- 
sults and others on symmetric monopoles will be presented elsewhere 0. 
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Note added. 

Recently the moduli space metric for the tetrahedrally symmetric 4-monopoles intro- 
duced in this paper has been calculated, P.M. Sutcliffe, Phys. Lett. 357B, 335 (1995). 
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